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Abstract. We present a new assessment of the contribution of the Blazar population to the extragalactic back- 
ground radiation across the electromagnetic spectrum. Our calculations rely on deep Blazar radio counts that we 
have derived combining several radio and multi-frequency surveys. We show that Blazar emission integrated over 
cosmic time gives rise to a considerable broad-band non-thermal Cosmic Background that in some parts of the 
electromagnetic spectrum dominates the extragalactic brightness. 

We confirm that Blazars are the main discrete contributors to the Cosmic Microwave Background where we 
estimate that their integrated emission causes an apparent temperature increase of 5-50 fj,K in the frequency range 
50-250 GHz. The CMB primordial fluctuation spectrum is contaminated starting at multipole I ~ 300 — 600, in the 
case of a completely random source distribution, or at lower I values if spatial clustering is present. We estimate 
that well over one hundred-thousand Blazars will produce a significant signal in the maps of the upcoming Planck 
CMB anisotropy mission. Because of a tight correlation between the microwave and the X-ray flux these sources 
are expected to be X-ray emitters with flux larger than a few 10 -15 erg cm~ 2 s~ 1 in the soft X-ray band. A large 
fraction of the foreground sources in current and near future CMB anisotropy maps could therefore be identified 
and removed using a mult i- frequency approach, provided that a sufficiently deep all sky X-ray survey will become 
available in the near future. 

We further show that Blazars are a major constituent of all high energy extragalactic backgrounds. Their con- 
tribution is expected to be 11-12 % at X-ray frequencies and possibly 100% in the ~ 0.5 — 50 MeV band. At 
higher energies (E > 100 MeV) the estimated Blazar collective emission, obtained extrapolating their integrated 
micro-wave flux to the 7-ray band using the SED of EGRET detected sources, over-predicts the extragalactic 
background by a large factor, thus implying that Blazars not only dominate the 7-ray sky but also that their 
average duty cycle at these frequencies must be rather low. Finally, we find that Blazars of the HBL type may 
produce a significant amount of flux at TeV energies. 
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1. Introduction 



The identification of the first quasar I Schmidt (1963) I 



marked the beginning of AGN astrophysics but also the 
discovery of the first flat spectrum radio quasar (3C273), 
a type of highly variable, often polarized extragalactic ra- 
dio sources that, together with the even more puzzling BL 
Lacertae objects, make the class of Blazars, the most ex- 
treme type of Active Galactic Nuclei known. The Blazar 
typical observational properties include the emission of 
electromagnetic radiation across the entire spectrum, from 
radio waves to the most energetic 7-rays, irregular rapid 
variability, apparent super-luminal motion, flat radio spec- 
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trum, large and variable polarization at radio and, espe- 
cially, at optical frequencies. 

Blazar of the BL Lacertae type (BL Lacs) are dis- 
tinguished by non-thermal emission with no (or very 
weak) emission lines, are often associated with the nu- 
clei of elliptical galaxies and are the only population 
of extragalactic sources that shows negative cosmolog- 
ical evolution (Bade et al.(1998)[ |Giommi et al. (1999)1 
Rector et al.(2000']^ Flat Spectrum Radio Quasars 
(FSRQs) share the strong non-thermal emission of BL 
Lacs but also show intense broad line emission and 
strong cosmological evolution similar to that of radio 
quiet QSOs (|Caccianiga et al. (2002)1 |Landt et al. (2001)1 
Wall fc Peacock(1985)|). " 
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Blazars are widely assumed to be powerful sources 
emitting a continuum of electromagnetic (e.m.) radiation 
from a relativistic jet viewed closely along the line of 
sight (|Blandford fc Rees(1978)| |Urry fc Padovani(1995) I. 



The broad-band electromagnetic spectrum is composed 
of a synchrotron low-energy component that peaks (in 
a Log{vf(v)) — Log{v) representation) between the far 
infrared and the X-ray band, followed by an Inverse 
Compton component that has its maximum in the hard 
X-ray band or at higher energies, depending on the lo- 
cation of the synchrotron peak, and extends into the 7- 
ray or even the TeV band. Those Blazars where the syn- 
chrotron peak is located at low energy are usually called 
Low energy peaked Blazars or LBL, while those where 
the synchrotron component reaches the X-ray band are 
called High energy peaked Blazars or HBL [see Fig. ^ and 
Padovani & Giommi(1995) . LBL sources are the large 
majority among Blazars [e.g., Padovani et al.(2003a)| and 
are usually discovered in radio surveys, while HBL objects 
are preferentially found in X-ray flux limited surveys since 
at these frequencies they are hundreds, or even thousands, 
of times brighter than LBLs of similar radio power. 
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Fig. 1. The Spectral Energy Distribution of differ- 
ent types of Blazars represented by Synchrotron-Self- 
Compton models with emission peaking at different en- 
ergies. LBL sources are those where the synchrotron peak 
is located at Infra-Red frequencies and therefore their X- 
ray emission is due to Inverse Compton radiation. HBL 
Blazars are those where the synchrotron component peaks 
in the UV/X-ray band making their / x //r flux ratio hun- 
dreds or thousands of times higher than that of LBLs. 

Radio quiet AGN are approximately one order 
of magnitude more abundant than Blazars and have 
been shown to be the major constituent of the 
Cosmic X-ray Back ground (CXB) (|Giacconi et al.(1962)| 
Rosati et al.(2002)[ |Moretti et al.(2003)| ) leading to the 
now consolidated picture in which the CXB is composed 
of radiation generated by the accretion onto super-massive 
black holes integrated over cosmic time. 



Despite the relatively low space density of Blazars, 
their strong emission across the entire electromagnetic 
spectrum makes them potential candidates as significant 
contributors to extragalactic Cosmic Backgrounds at fre- 
quencies where the accretion mechanism does not produce 
much radiation. These extragalactic backgrounds would 
then be mostly composed of non-thermal radiation gener- 
ated in Synchrotron/Inverse Compton-type environments. 



Recently, Giommi & Colafrancesco(2004) showed that 
Blazars are by far the largest population of extra- 
galactic objects detected as foreground sources in CMB 
anisotropy maps and that their emission contaminates 
the CMB angular power spectrum at a significant level. 
Padovani et al.(1993)[ based on the detection of a small 
sample of Blazars at 7-ray frequencies concluded that the 
Blazar population should produce a large fraction of the 
high energy Cosmic Background. 

In this paper we re-assess the Blazar contribution to 
the Cosmic energy in the microwave (CMB), in the X- 
ray (CXB), in the Gamma-ray (CGB) and TeV (CTB) 
part of the electromagnetic spectrum. Our calculations 
rely on a new estimation of the Blazar radio LogN- 
LogS, that we have assembled combining several radio 
and multi-frequency surveys, on flux ratios in different en- 
ergy bands and on observed Blazar broad-band Spectral 
Energy Distributions (SED). 

2. The Blazar radio LogN-LogS 

In this section we use several radio and multi-frequency 
surveys to build a new deep Blazar radio LogN-LogS that 
updates and extends to lower fluxes the counts presented 
in Giommi & Colafrancesco(2004) Since the surveys con- 



sidered have been carried out at three different observing 
frequencies (1.4, 2.7 and 5 GHz), we convert all flux den- 
sities to a common frequency before proceeding. We take 
5 GHz as the reference frequency and (unless otherwise 
stated) we apply the flux conversions assuming a spectral 
slope cv r = 0.25 (/ ex v~ ai ) which is approximately equal 
to the average value in all our samples. 

The derived Blazar counts are shown in Fig. [21 and 
can be described by a broken power law with parameters 
defined in the following equation: 



N(> S) = 



5.95 1CT 3 x S- 
0.125 x S~°- 9 



1.62 



S > 0.015 Jy 
S < 0.015 Jy 



(1) 



The slope up to 15 mJy (dot-dashed line) is a good rep- 
resentation of the data, whereas the flattening at fainter 
fluxes (dotted line) is necessary to avoid that the pre- 
dicted Blazar space density exceeds the total density of 
NVSS radio sources at a few mJy. The slope below the 
break is somewhat arbitrary as only upper limits are avail- 
able in this flux regime. We have chosen a value of 0.9 
since this is the average slope of the LogN-LogS of ra- 
dio quiet AGN in the two flux decades below the break 
( Rosati et al.(2002)| |Moretti et al.(2003)| ) and is consis- 
tent with all the available constraints. 
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In the following we discuss the details of all the surveys 
in order of decreasing radio flux limit. 

2.1. Two Jansky flat spectrum radio survey 



The two Jansky 2.7 GHz sample (Wall fc Peacock(1985) I 
is based on a complete radio flux limited survey of 
flat spectrum (a r < 0.5, / oc ^ _Qr ) sources cover- 
ing the entire sky with the exclusion of the Galactic 
plane (|6| > 10). The sample includes 60 Blazars 



( di Serego-Alighieri et al.(1994)[|Urry fc Padovani(1995) \ 
corresponding to a space density of 0.002 deg~ z . This 
value has been plotted as an open square symbol in Fig. 
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2Jy Hot Spectrum Radio Survey 

Uy ASDC-RASS-NVSS 

WMAP selected Blazars 
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Fig. 2. The radio (5GHz) LogN-LogS of Blazars built 
combining several radio and multi-frequency surveys. The 
different symbols represent the surveys used (see text for 
details). The filled square symbols in the lower part of 
the diagram represent the LogN-LogS of extreme HBL 
BL Lacs from the Sedentary Survey. Although these ob- 
jects are a tiny minority (« 0.1%) of the overall Blazar 
population they play an important role in producing the 
Cosmic Background at X-ray and at very high energies 
(7-ray/TeV). 



2.2. One Jansky ASDC-RASS-NVSS Blazar Sample 

The ASDC-RASS-NV SS Uy Blazar (Uy-ARN) Survey 
( Giommi et al. (2002c I is a radio flux limited (f r > 
lJy@1.4GHz ) sample of Blazars selected by means of 
a cross correlation between the ROSAT All Sky Survey 
(RASS) catalog of X-ray sources ( |Voges et al.(1999)|) and 
the subsample of NVSS survey ( Condon et al. (1998) ) 
sources with flux density larger than 1 Jansky. The 
selection was carried out applying the same multi- 
frequency technique used for the definition of the ex- 
treme HBL BL Lac sample in the Sedentary survey 
(Giommi et al.(1999)|). To avoid the complications due to 



plane, only sources located at Galactic latitudes larger 
than 20 degrees were considered. The accurate radio po- 
sitions allowed us to obtain firm associations with all the 
optical counterparts. The sample includes a total of 226 
sources, 96% of which were previously known objects, in 
most cases well documented in the literature as they are 
bright radio sources. The class composition of the sample 
is reported in Table 1. 



Table 1. One Jansky ASDC-RASS-NVSS sources 



Number 
of sources 



Classification 



132 FSRQs 

25 BL Lac 

36 SSRQs 

21 radio galaxies 

3 other galaxies 

9 unidentified (3 flat, 6 steep) 



From past X-ray measurements we know that 
Blazars are characterized by a / x // r ratio ranging 



from 
10" 9 



1. x 10 



-13 



erg cm s" 



to over 1. 



high Nh photoelectric absorption close to the Galactic 



erg cm _2 s _1 Jy _1 ( Padovani(2002) \ and therefore 
all Blazars with flux larger than 1 Jy should be de- 
tectable in a survey with X-ray sensitivity of 1. x 10 -13 
erg cm~ 2 s _1 Jy -1 . However, the X-ray limiting flux of the 
RASS strongly varies across the sky depending on the ef- 
fective exposure and on the amount of Galactic absorption 
(Nh) along the line of sight. 

Figure shows the high Galactic Latitude (\b\ > 20) 
sky coverage of the RASS in the region of overlap with 
the NVSS survey. The curve was calculated taking into 
account of the ROSAT effective exposure in sky bins of 1 
square degrees in size and converted to sensitivity assum- 
ing a power law spectrum with energy index a — 1.0 and 
setting the amount of Nh equal to the Galactic value in 
the direction of each sky bin. 

The X-ray sensitivity ranges from somewhat below 1 x 
10 -13 erg cm _2 s~ 1 , where about 1000 square degrees of 
sky are covered, to « 2. x 10 -12 erg cm _2 s _1 , where the 
maximum area of 23,000 square degrees is reached. 

Assuming that the three unidentified flat spectrum 
sources are Blazars, the number of objects of this type 
(FSRQ+BL Lacs+unidentified sources with cv r < 0.5) in 
the Uy-ARN survey is 160. 

This sample has been used to estimate the Blazar space 
density above Uy taking into account that the sky cover- 
age of Fig. implies that faint sources with flux around 
e.g. 1. x 10 -13 erg cm~ 2 s _1 could be detected only in 
« 1,000 square degrees of sky, whereas brighter X-ray 
sources with flux of e.g. 1. x 10~ 12 erg cm _2 s _1 could be 
detected over a much larger portion of the sky (~ 20, 000 
square degrees). The counts have been converted to 5GHz 
and then plotted as open circles in Fig. [3 
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Fig. 3. The high Galactic Latitudes (\b\ > 20) sky cover- 
age of the RASS survey in the region of overlap with the 
NVSS catalog (Dec > -40 degrees). 



2.3. WMAP microwave selected Blazars 

The Wilkinson Microwave Anisotropy Probe (WMAP) is 
a space observatory dedicated to the accurate investiga- 
tion of primordial fluctuations in the Cosmic Microwave 
Background (Bennett et al. (2003a)). A catalog of 208 



bright foreground sources detected in one or more of the 
five microwave WMAP channels during the first year all 



sky survey has been published by Bennett et al. (2003b) 
With very few exceptions all entries are well known bright 
sources at cm wavelengths and in most cases observed at 
several radio frequencies. Giommi & Colafrancesco(2004) 
have recently shown that the large majority of these ob- 
jects are Blazars. 

To determine how WMAP selected Blazars contribute 
to the 5GHz counts we have considered all high Galactic 
latitude (|6| > 20 deg) detections with flat radio spectrum 
(cv r < 0.5, / oc v~ ai ) and with SED typical of Blazars. 
Table |2 gives the detailed statistics: 

Table 2. WMAP bright foreground source catalog 



Number 


Classification 


of sources 




141 


FSRQs 


23 


BL Lac 


13 


radio galaxies 


5 


SSRQs 


2 


starburst galaxies 


2 


planetary nebulae 


17 


unidentified 


5 


no radio counterpart 




(probably spurious) 



The space density of WMAP detected sources with 
radio flux above 1, 3, 5 and 10 Jansky at 5GHz has been 
plotted in Fig. [3 as filled star symbols. 



With the exception of the point at 1 Jansky, 
which is most likely an underestimation due to in- 
completeness of the WMAP catalog at this flux limit 
(Bennett et al. (2003b) I, the agreement with other radio 



surveys at cm wavelengths is very good implying that the 
Blazar selected by WMAP and at radio frequencies belong 
to the same population of objects. 

2.4. The Parkes quarter-Jansky Flat Spectrum sample 

The Parkes quarter-Jansky Flat Spectrum sample 



( Wall et al.(2004)| ) is a 100% identified radio flux limited 
survey at a frequency of 2.7 GHz. 

A total of 328 FSRQ and 43 BL Lacs have been de- 
tected with flux density higher than 0.25 Jy in an area 
of 8785 square degrees (sample 1 of Wallet al.(2004)). 
Considering that this survey only accepts sources with 
spectral index flatter than a r < 0.4 (/ oc v~ ar ) we rescale 
the sample density by a factor 1/0.75 which is the ratio of 
QSOs and BL Lacs with a T < 0.7 to those with cv r < 0.4 
in the Parkes survey known as PKSCAT90 with flux den- 
sity larger than 0.25 Jy. A similar ratio is present in the 
Uy-ARN survey. 

The space density of Blazar in this survey is therefore 
0.06 objects per square degree. We convert 2.7GHz fluxes 
to 5GHz and plot the density as an open diamond in Fig. 

m 

We note that 84 additional sources are classified by 



Wall et al.(2004) as flat radio spectrum Galaxies (57 of 
which without redshift); as some of these may well be BL 
Lacs the Blazar content of this sample will probably grow 
in the future. 



2.5. The DXRBS Blazar Survey 

The Deep X-Ray Radio Blazar Survey (DXRBS) is a 
radio flux limited sample based on a double selection 
technique at radio and X-ray frequencies and uses op- 
tical data to refine the sample. DXRBS searches for 
Blazar candidates among serendipitous X-ray sources of 
the ROSAT PSPC pointed observations listed in the 



WGA catalog (White et al.(1995)l, restricting the sam- 
ple to objects with radio spectral index flatter than a r = 
0.5 and with broad-band spectral indices a TO and a ox in 
the region occupied by Blazars. Details of the selection 
method and of the optical identification of the candidates 
are described in Permian et al.(1998)[ Landt et al.(2001){ 



Padovani et al.(2005) 



Although the radio limit of the survey is 50 mJy at 
5GHz, the X-ray sensitivity of the ROSAT PSPC point- 
ings is not deep enough to ensure completeness at all radio 
fluxes, as some Blazars with flux lower than 100 mJy are 
expected to have an X-ray flux below the ROSAT WGA 
limit of ps 2 x 10~ 14 erg cm _2 s _1 . We have therefore es- 
timated the Blazar space density taking into account the 
latest LogN-LogS results of |Padovani et al.(2005) for both 
FSRQs and BL Lacs and correcting the points at 50 and 
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100 mJy for the fraction of lost objects as predicted from 
the distribution of Fig. 0] . 

The final Blazar space density from this survey is plot- 
ted in Fig. [2] as open triangles. 

2.6. The ASDC-XMM-Newton-NVSS (AXN) Sample 

In order to push the radio limit to fluxes significantly be- 
low 50 mJy at 5GHz using the same multi-frequency se- 
lection method of the DXRBS survey it is necessary to 
reach X-ray sensitivities proportionately deeper than that 
of the ROSAT-PSPC. For this purpose, we have searched 
for serendipitous NVSS radio sources in XMM-Newton 
EPIC-pn ( jStriider et al.(2001) I X-ray images which pro- 
vide approximately one order of magnitude better sensi- 
tivity than ROSAT. 

We have used the data in the XMM-Newton public 
archive available at the ASI Science Data Center (ASDC) 
in December 2004. At that date a standard processing, 
which removes periods of high background resulting from 
solar flares and detects serendipitous sources, had been 
run at ASDC on 1220 EPIC-pn fields, 847 of which at 
high galactic latitude (\b\ > 20). 

Of these EPIC-pn observations we have considered the 
subsample of 188 non-overlapping fields that satisfy the 
conditions listed below that are necessary to ensure that 
the sample of serendipitous sources is suitable for our pur- 
poses: 

1. the EPIC-pn instrument was used in full window 
mode; 

2. the pointing of the XMM observation was at declina- 
tion north of -40 degreed to ensure overlap with NVSS 
survey; 

3. no bright or extended source is present in the field of 
view; 

4. whenever more than one exposure of the same field was 
available in the archive, the deeper one was taken. 

In addition, we have excluded a 5 arc-minutes circu- 
lar area around the target of each field. The total area 
covered by the 188 fields is 26.3 square degrees. A careful 
inspection of all the X-ray/radio/optical associations with 
broad-band spectral indices cv ro and a ox in the region oc- 



cupied by Blazars I Giommi et al.(1999) I led to the results 
summarized in Table |3 



Table 3. NVSS radio sources detected in XMM-Newton 
EPIC-pn X-ray images 



Number 
of sources 


Space density 
sources deg" 1 


Flux limit 
1.4 GHz 


Flux limit 
5 GHz * 


50 


1.9 ±0.3 


50 


31 


107 


> 4.0 


20 


12 


149 


> 5.6 


10 


6 



For this first (conservative) estimation we assume that 
all XMM fields are sufficiently deep to detect all Blazars 
above 50 mJy. Because of the spread in exposure times and 
of the reduction in sensitivity at large off-axis angles due 
to vignetting effects and PSF degradation, some objects 
(of the order of 20-30%) might have been missed and the 
number of Blazars found is only a lower limit. On the other 
hand, from the Uy-ARN and other surveys, we expect 
that a similar percentage of candidates be associated to 
steep spectrum objects, so that the two effects roughly 
cancel out. We will discuss the impact of this reduction 



in sensitivity in a future paper ( Cavazzuti et al.(2005) I 
where the sky coverage of the sample will be fully taken 
into account. 

The space density at 50 mJy at 1.4GHz (31 mJy at 
5. GHz) has been plotted in Fig. [3| as a filled circle. The 
points at 12 and 6 mJy are drawn as lower limits because 
in this case a significant fraction (e.g. ^ 50%) of these 
faint Blazars (especially at 6 mJy) is expected to be below 
the sensitivity limit of our XMM-Newton X-ray images. 
In fact, the / x //r distribution of known Blazars (see Fig. 

implies that the expected soft X-ray flux of a 10 mJy 
Blazar is below 2 x 10~ 15 erg cm" 2 s _1 in about 30% of 
the cases. These X-ray sources could not be detected in 
most of our EPIC images, especially at large off-axis angles 
where most of the area is located. 

2.7. The Sedentary Survey of extreme High Energy 
Peaked BL Lacs 

We conclude this section summarizing some of the 
results from the multi- frequency 'Sedentary Survey' 
( Giommi et al.(1999)[ 



Piranomonte et al.(2004) 



We assume a spectral slope a r =0.4 



Giommi et al. (2005) I which is a deep (f r > 3.5 mJy 
at 1.4 GHz), 100% identified radio flux limited sample of 
150 extreme HBL objects characterized by a / x //r ratio 
higher than 3 x 10~ 10 erg cm" 2 s" 1 Jy" 1 . This survey does 
not have a direct impact on the full Blazar LogN-LogS 
of Fig. |21 since such extreme HBL objects only represent 
a tiny minority of the full Blazar population at radio 
frequencies. However, the very high / x //r flux ratio, 
which in these sources ranges from five hundred to over 
five thousand times that of typical Blazars, makes these 
rare radio emitters a numerically important population 
of sources at X-ray and TeV frequencies and therefore 
makes them potentially significant contributors to the 
Cosmic Background radiation at these frequencies. The 
space density of the objects at fluxes f r > 6 mJy (5 GHz) 
in this survey is reported in Fig. [21 as filled squares. 



3. From radio frequencies to other spectral bands 

Once the LogN-LogS of a population of sources is known 
in a given energy band, it is possible to estimate their 
emission in other parts of the electromagnetic spectrum 
provided that the flux ratio in the two bands, or even 
better, the overall energy distribution, is known. 
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In this section we deal with flux ratios and SEDs of 
Blazars that will be used later in the paper to estimate 
the contribution to Cosmic Backgrounds at frequencies 
higher than the radio ones. 

The collection of SED of a large sample of Blazars built 
using BeppoSAX and multi-frequency data presented by 



Giommi et al. (2002a) clearly shows that a wide variety of 
broad-band spectral forms exist. This large spectral diver- 
sity, however, can be reproduced, at least in first approx- 
imation, by Synchrotron Self Compton (SSC) emission 
models as those shown in Fig. ^ Sometimes, to explain 
the 7-ray part of the distribution, it is necessary to add 
other components to this simple picture. We will deal with 
this possibility later when we consider the contribution of 
Blazars to the 7-ray background. 

Although a continuous range of synchrotron peak en- 
ergies probably exists, it is useful to divide the population 
of Blazars into Low energy peaked Blazars (or LBL) whose 
X-ray emission is due to Inverse Compton radiation, and 
in High energy peaked Blazars (HBL) characterized by 
a much higher X-ray flux (for the same radio flux) since 
their synchrotron component extends into the X-ray band. 
This simple separation in two broad categories also reflects 
the early selection methods since LBL objects have mostly 
been discovered in radio surveys, while HBL objects have 
been typically found in X-ray flux limited samples. 

The large difference in X-ray emission in LBL and 
HBL objects is graphically reflected in the distribution 
of f x / f T flux ratios shown in Fig. 0] which spans about 4 
orders of magnitudes. 



0.3 



0.25 



3 0.15 




0.05 



Fig. 4. The / x // r distribution of Blazars estimated from 
the Uy-ANR and from a sample of about 2000 HBL 
Blazar candidates (see text for details). 



The plot shown in Fig. 0] was built using two ra- 
dio flux limited samples. At / x // r values lower than 
5 x 10 -12 erg cm _2 s _1 Jy _1 (corresponding to LBL ob- 
jects) we have used the Uy-ARN sample corrected for 
the RASS sky coverage. As the percentage of Blazars 
sharply drops at / x // r values larger than 5 x 10 -12 
erg cm _2 s _1 Jy _1 the Uy-ARN sample rapidly becomes 
statistically inadequate and lager samples are clearly 
needed. We have therefore built a radio flux limited sam- 
ple including nearly 2000 Blazar candidates with / x / / r > 
5 x 10 -12 erg cm~ 2 s~ 1 Jy~ 1 applying the same technique 
used for the Sedentary survey which is about 85% effi- 
cient in selecting HBL Blazars with / x // r > 3 x 10 -10 
erg cm _2 s _1 Jy _1 ( jGiommi et al.(2005) l. 



To that purpose, we have cross-correlated the posi- 
tions of the X-ray sources in the ROSAT All Sky Survey 
( |Voges et al.(1999)|l with the radio sources of the NVSS 
catalog (Condon et al. (1998) I and we have estimated op- 
tical (J or F) magnitudes from the GSC2 Guide Star 
Catalog (McLean et al.(2000)l. We have then calculated 
the broad-band a ox and a ro spectral slopes and we have 
only accepted objects in the area of the a OK - a TO plane 



occupied by Blazars (Giommi et al.(1999) I. Whenever no 
optical counterpart was found within the radio positional 
uncertainty we have assumed a lower limit of Jmag = 19.5, 
which is the limit of the GSC2 catalog. 

To estimate the distribution at / x // r > 5 x 10~ 12 
erg cm~ 2 s _1 Jy _1 the radio flux limit of the sample must 
correspond to an X-ray flux at which the ROSAT sur- 
vey covers a sky area large enough to allow the detection 
of at least a few objects. We have chosen f r =25 mJy, 
corresponding to fx > 1-25 x 10 -13 erg cm _2 s _1 where 
the RASS covers about 2000 square degrees of sky. We 
have then derived the tail of the / x // r distribution for 
/ x //r < 5 x 10~ 12 erg cm _2 s _1 Jy~ 1 taking into account 
the RASS sky coverage shown in Fig. [21 



We have checked the reliability of this selec- 
tion method using the subsample of 514 objects 
for which Sloan Digital Sky Survey (SDSS) [e.g., 
York et al. (2000)1 |Stoughton et al.(2002)| spectral data 



( Turriziani et a l. (2005) ) are available and we have found 
that over 80% of the 319 candidates with / x // r > 5 x 10~ 12 
erg cm _2 s _1 Jy _1 are indeed spectroscopically confirmed 
Blazars. 

To properly normalize this distribution we scaled it 
so that the fraction of Blazars with / x // r > 3 x 10 -10 
erg cm~ 2 s -1 Jy -1 (the / x // r limit of the Sedentary sur- 
vey) is equal to the density ratio between the full Blazar 
population and that of the extreme HBL of the sedentary 
survey (dashed line and filled squares in Fig. El- 

The combined / x // r distribution is shown in Fig. 0] 
where we see that the large majority (94%) of Blazars are 
of the LBL type (defined here as objects with / x // r < 
5 x 10~ 12 erg cm _2 s _1 Jy _1 ). 
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Fig. 5. The extragalactic Cosmic Background Energy dis- 
tribution at microwave, X-ray and 7-ray energies. 



4. The contribution of Blazar to Cosmic 
Backgrounds 

Figure [S] shows the spectral energy distribution of the 
Cosmic Extragalactic Background radiation in the mi- 
crowave, X-ray and 7-ray band, that is where we expect 
that the Blazar collective emission gives a significant con- 
tribution. 

The Cosmic Microwave Background (CMB), taken 
from the very accurate measurements of the COBE satel- 
lite is represented by a black-body spectrum with tem- 
perature of 2.725 °K (|Mather et al.(1999)|>; the X-ray 



Background (CXB) is taken from HEAO-1 measure- 
ments ( |Marshall et al.(1980)| |Gruber et al.(1999)| > and 
has been scaled to match the more recent BeppoSAX , 
ASCA and X MM-Newton results in the softer 2- 
10 keV band (|Vecchi et al.(1999)1 |Lumb et al.(2002)| 
Kushino et al.(2002)l. The Gamma Ray Background is 
derived from the COMPTEL data in the range 0.8 - 30 



MeV < |Kappadath(1998)| ), and from EGRET data in the 
range 30 MeV - 50 GeV. We recall that the estimate of 
the extragalactic 7-ray background emission depends on 
the Galactic diffuse emission model which itself is not 
yet firmly established. In Fig. we report the results 
of two different analysis of the EGRET data: open cir- 
cles from |Sreekumar et al.(1998) and filled circles from 
Strong et al.(2004)[ which uses an improved model of the 
Galactic diffuse continuum gamma-rays. As for the TeV 
diffuse background, we report the upper limits in the 20 - 
100 TeV region derived from the HEGRA air shower data 
analysis ( |Aharonian et al. (2002)1 ). Note that the HEGRA 
measurements are sensitive both to non-isotropic (galac- 
tic) and isotropic (extragalactic) component. In the 1 TeV 
- 1 PeV energy range, other experiments give only upper 
limits and there is no clear observation of a diffuse photon 
signal yet. 

In the following we estimate the Blazar contribution 
to the above described Cosmic Backgrounds basing our 



calculation on the radio LogN-LogS built in Sect. 2 and on 
flux ratios in different bands or on the observed Spectral 
Energy Distributions of 7-ray detected sources. 

4.1. The Cosmic Microwave Background 

The contribution of Blazars to the Cosmic Microwave 
Background (CMB) has been estimated in the past 
from different viewpoints (see, e.g., Toffolatti et al.(1998) 
Giommi & Colafrancesco(2004) I. 

Here we use the Blazar radio LogN-LogS 
of Fig. El to update the results presented in 



Giommi & Colafrancesco(2004) 

From Fig. 01 we see that the Blazar LogN-LogS can be 
represented by a broken power law model with alpha=1.62 
(integral slope) up to a break point where the slope flat- 
tens significantly. The precise position of the break and 
the amount of flattening cannot be estimated with our 
data. However, a break must occur somewhere around 10- 
15 mJy otherwise the number of Blazars would exceed the 
total number of radio sources at about 2-3 mJy. In the fol- 
lowing we use the somewhat conservative values of 15 mJy 
for the break and 0.9 for the slope below the break. Both 
values are consistent with the lower limits on the Blazar 
counts at 12 and 6 mJy (see Sect. 2) 

The integrated background intensity due to Blazars 
therefore can be expressed as 



IbIciz 



Uy 
0.1m Jy 



Q dN rlQ 



(2) 



where dN/dS is the differential of eq. 0). The minimum 
integration flux of 0.1 mJy for S m i n is likely to be conser- 
vative since Blazars with radio flux near or below 1 mJy 
are already included in the Einstein Medium Sensitivity 
Survey BL Lac sample (Rector et al. (2000) 1. The inte- 
grated intensity Isiazars is then converted from 5GHz 
to microwave frequencies convolving the flux value with 
the observed distribution of spectral slopes between 5GHz 
to microwave frequencies as estimated from the Uy-ARN 
sample (se eJGiommi fc Colafrancesco( 200*4) for details). 

Figure El plots the fractional contamination, defined 
as iBiazars/ Icmb, that is the ratio between the Blazar 
integrated emission and the CMB intensity (solid line), 
and the equivalent apparent temperature increase of the 
CMB (dotted line), as a function of frequency. 

The contribution of Blazars to the temperature 
anisotropy spectrum of the CMB is calculated as 
(AT) Blazar = [(2n)- 1 C e £(e+ l)] 1 / 2 , where 



Ct 



Blazar 



(3) 



assuming 
tributed 



that 
like 



(Tegmark fc Efstathiou(1996) 



the Blazars are spatially dis- 
a Poissonian random sample 
Scott fe White(1999)| 



Giommi & Colafrancesco(2004) 1. The quantity at the 
right hand side in eq. (j2J) is the usual Poisson shot-noise 
term ( |Peebles(1980)l |Tegmark fc Efstathiou(1996)l ). We 
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Fig. 6. The Blazar fractional contamination of the CMB 
(solid line) and the apparent increase in CMB temperature 
due to Blazars (dotted line) as a function of frequency 

neglect here the clustering term, w(/) 2 , which adds 
to the Poissonian term since there is not yet a clear 
estimate of the Blazar clustering. We note however, that 
the inclusion of clustering can significantly increase the 
amount of contamination, especially at large angular 



scales ( Gonzales-Nuevo et al.(2005)i. 
For the Blazar population described by the LogN-LogS 
given in fig. we found Blazar ~ 54.4 Jy 2 sr _1 at 
41 GHz and Q >B i azar « 51.6 Jy 2 sr" 1 at 94 GHz. Our 
results can be translated into temperature units using the 
conversion between the isotropic black-body (Planckian) 
brightness Bq{v) and the CMB temperature To which 
can be written as 



dB 
dT 



kTo 
he 

24.8Jy sr" 



sinh(x/2) 



sinh(x/2) 



(4) 



Here x = hv/kT — ^/56.84 GHz is the a-dimensional 
frequency given in terms of the Planck constant h, of the 
speed of light c and of the Boltzmann constant k. We 
found a value C^Biazar ~ 2.22 • 10~ 2 /xK 2 sr at 41 GHz 
and C^Biazar ~ 1-09 • 10~ 3 yK 2 sr at 94 GHz. We show 
in Figs. and |H1 the quantity ATeiazar and we compare 
it to the CMB fluctuation power spectrum which best 
fits the available data. The previous value is shown as 
solid lines and has to be considered as a definite lower 
limit for C^Biazar since it neglects the contribution of 
steep-spectrum sources at low frequencies which flatten 
at these frequencies (41 and 94 GHz) and the effect of 
flux variability. The dashed lines in Figs. [7| and [S] show 
the Blazar power spectrum by adding the contribution 
of radio sources with steep-spectrum at low radio fre- 
quencies which flatten at higher frequencies. The effect 
of spectral and flux variability also allows for an increase 
of the fluctuation level C^Biazar because many Blazars 




Fig. 7. The contribution of Blazars to the CMB fluctu- 
ation spectrum in the WMAP Q channel at 41 GHz as 
evaluated from the LogN-LogS given in Fig. [3 (solid line). 
We also show the angular power spectrum for the Blazar 
population by adding an estimate of the possible contribu- 
tion of radio sources with steep-spectrum at low radio fre- 
quencies which flatten at higher frequencies (dashed line) . 
The dotted line includes also the effect of spectral and 
flux variability (see text for details). Although this addi- 
tional contamination may be substantial a precise estima- 
tion can only be done through simultaneous high resolu- 
tion observations at the same frequency. A typical CMB 
power spectrum evaluated in a ACDM cosmology with 
Q m = 0.3, 51a = 0.65, fib = 0.05 which best fits the avail- 
able data is shown for comparison. 



below the sensitivity threshold of CMB experiments can 
be detected during flares. The Blazar flux variability at 
millimeter wavelengths may be very substantial (higher 
than factors 3-10 on time scales of weeks to years seen at 
cm wavelengths) and could definitely increase the level of 
contamination of CMB maps when these are built over 
long integration periods. The dotted lines in Figs. and 
|H1 show the Blazar power spectrum when also this last 
contribution is taken into account. 

Our numerical estimates are subject to some uncer- 
tainties which depend on the parameters of the LogN- 
LogS that are not well constrained. For instance if instead 
of using the conservative values of 15 mJy and 0.9 for the 
break flux and slope below the break we used the less con- 
servative values of 10 mJy and a = 1.0 we would have an 
increase of « 50% for I B iazars and of «4% for CYsiazar • 

The contamination level shown in Figs.[7|and|Hlis that 
expected in case no Blazars are removed from the CMB 
data. The observed WMAP angular power spectrum re- 
ported by Hinshaw (2003) has been obtained using a pro- 
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the range 20-25 %. This estimate does not include possi- 
ble variability effects and additional core-dominated radio 
sources 




WMAP 
94 GHz 



100 
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Fig. 8. The contribution of Blazars to the CMB fluctua- 
tion spectrum as in Fig.0for the WMAP 94 GHz channel. 



cedure in which a pure primordial signal has been initially 
assumed, and then subtracting the foreground sources 
that could be recognized as such in the WMAP data, 
thus removing only the bright tail of the Blazar LogN- 
LogS distribution. The procedure adopted to derive the 
CMB power spectrum should take into account the full 
point-like source contribution implied by our LogN-LogS. 
Such an approach would both influence the shape of the 
expected power spectrum and increase the statistical un- 
certainties of the WMAP data points, especially at high 
multipoles, where the Blazar contribution is larger. 

We stress that the previous calculations are performed 
neglecting the clustering term oj(I) 2 and thus they must 
be considered again as a lower limit to the realistic angu- 
lar power spectrum contributed by the Blazars. The effects 
of clustering of FSRQ on the CMB fluctuation spectrum 
has been partially estimated by some authors: for instance 

calculated from simulations 



Gonzales- Nuevo et al.(2005) 



that the clustering of extragalactic radio sources is rather 
small at i = 200 for the Planck frequencies, while the con- 
tribution of the clustering term to the confusion noise is 
likely to be the dominant one. Scott and White (1999) es- 
timated instead that the clustering of SCUBA sources to 
the expected Planck CMB spectrum provides an increase 
of a factor 5-10 with respect to the Poissonian term for 
the point-like source contribution at I = 200. It is clear 
that the expectations for the clustering effect strongly de- 
pend on the adopted model for the source counts and on 
their clustering model. Based on the correlation function 
of Loan et al.(1997) (adopted by Toffolatti) and on that 
for the SCUBA sources adopted by |Scott h White(1999j| 
and using our Blazar LogN-LogS, we expect that the con- 
tamination of the first peak of the fluctuation spectrum 
(at the WMAP 41 GHz channel) is at a level comprised in 



4.2. The Soft X-ray Background 

We estimate the contribution of Blazars to the CXB at 1 
keV using two methods: i) converting the integrated ra- 
dio flux calculated with eq. into X-ray flux using the 
observed distribution of / x // r flux ratios of Fig.0J and ii) 
converting the integrated Blazar contribution to the CMB 
(at 94 GHz) to X-ray flux using the distribution of a^x , 
the microwave (94 GHz) to X-ray (IkeV) spectral slope, 
defined as 



Logjflkev/ fviGHz) _ Log(fikev/ fviGHz) 

Log(vi keV /is 9iGHz ) 



6.41 



(5) 



and estimated from the subsample of Blazars detected by 
WMAP in the 94 GHz channel and for which an X-ray 
measurement is available (see Fig.|UJ). 

The first method gives a total Blazar contri- 
bution to the X-ray background of 2.7 x 10~ 12 
erg cm _2 s _1 deg _1 (about 70% of which due to HBL 
sources with / x // r > 5 x 10 -12 erg cm _2 s _1 Jy _1 ) in the 
ROSAT 0.1-2.4 keV energy band. Assuming an average 
Blazar X-ray energy spectral index of a x — 0.7 (which is 
a compromise between a flatter spectrum that is typical 
of LBL sources and an cv x steeper than 1 for HBL objects) 
this flux converts to 2.6 x 10 -12 erg cm _2 s _1 deg _1 in the 
2-10 keV band or 11% of the CXB which is estimated to be 



2.3 x 10~ n erg cm _2 s -1 deg _1 (Perri fc Giommi(2000) ) 




Fig. 9. The distribution of spectral slopes between the mi- 
crowave (94 GHz) and X-ray (IkeV) band measured from 
the sample of 42 Blazars included in the WMAP catalog of 
bright microwave sources detected in the 94 GHz channel 
and for which X-ray measurements are available. 

The distribution of microwave to X-ray spectral slope 
of Fig. [3] has an average value of ( a^ x ) — 1.07 and a 
standard deviation of 0.08 corresponding only to about a 
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factor 3 in flux. This distribution is much narrower than 
that between the radio and X-ray band (Fig. 0] / x // r ) 
and the dispersion is comparable to the variability that 
Blazars show at radio and especially at X-ray frequencies. 

Since the 94 GHz and X-ray measurements are not 
simultaneous most of the dispersion in the distribution 
shown in Fig. [5] is probably due to variability both in 
the microwave and X-ray band. The intrinsic dispersion is 
therefore likely to be much smaller. For this reason in con- 
verting from the background level at 94 GHz estimated in 
the previous section to 1 keV we use the average a^ x value 
and assume no dispersion. We must note however that, 
just like the Uy radio surveys, the WMAP a^ x distribu- 
tion hardly includes any HBL sources, which make only a 
few percent of the population. 
From eq. JSJ) we get 

jW = 1-4 x 1CT 7 f 94GHz ■ (6) 

Since the Blazar integrated emission at 94 GHz 
is 7.2 x 10- 6 CMB 94G _ f/z or 0.64 Jy/deg 2 and 
the Cosmic X-ray Background is 2.3 x 10~ n 
erg cm^s^ 1 ( |Perri fc Giommi(2000)| ) in the 2-10 
keV band, equivalent to 2.31/xJy/deg 2 at 1 keV, eq. © 
for fg4GHz=0M Jy/deg 2 gives h ke v = 0.09 /iJy/deg 2 or 
3.9% of the CXB. 

Taking into account that the distribution of Fig.[§|only 
includes LBL objects and that HBL sources make two 
thirds the total contribution to the CXB (see above), the 
percentage scales to about 12% which is very close to the 
11% obtained with the previous method. 

Both result are in good agreement with the indepen- 



dent estimate of Galbiati et al.(2004) based on the XMM- 
Newton Bright Serendipitous Survey, who conclude that 
the radio loud AGN content of the CXB is 13%. 

4.3. Hard X-ray - soft j-ray Background 

The number of sources detected at energies higher than 
soft X-rays is still rather low and building reliable distri- 
butions of flux ratios between radio or microwaves and the 
Hard X-ray/7-ray fluxes similar to those of Fig. 0] and El is 
not currently possible. In order to estimate the Blazar con- 
tribution to high energy Cosmic Backgrounds (E > 100 
keV) we have therefore followed a different approach. We 
have extrapolated the predicted Blazar integrated inten- 
sity at microwave frequencies (eq.|2J) to the hard X-ray and 
soft 7-ray band using a set of hypothetical Synchrotron 
Sclf-Compton Spectral Energy Distributions. 

Figure HO shows the CMB, CXB and CGB together 
with three predicted SED from a simple homogeneous SSC 
models whose parameters are constrained to 1) be con- 
sistent with the expected integrated flux at 94 GHz, 2) 
have the a^ x slope equal to the mean value of the WMAP 
Blazars (a^ x — 1.07) and 3) possess a radio spectral slope 
equal to the average value in the WMAP sample. The 
three curves, forced to pass through the three star sym- 
bols graphically representing the three constraints listed 
above, are characterized by synchrotron peak frequencies 



of v peak = 10 12 - 8 , 10 13 - 5 and 10 13 8 Hz. From Fig.UHIwe see 
that a high value of v pea k over-predicts by a large factor 
the observed Hard- X-ray to soft 7-ray (w 10 20 — 2 x 10 22 
Hz or k, 500 keV-10 MeV) Cosmic Background, whereas 
a too low value of v pea k predicts a negligible contribution. 
The case v pea k = 10 13 ' 5 Hz predicts 100% of the Cosmic 
Background. Since the Log(^ pea fe ) values of Blazars in 
the Uy-ARN survey and WMAP catalog peak near 13.5 
and range from 12.8-13.7 within one sigma from the mean 
value, we argue that the data presently available indi- 
cate that Blazars may be responsible for a large frac- 
tion, possibly 100%, of the Hard-Xray/Soft 7-ray Cosmic 
Background. 
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Fig. 10. The possible contribution of LBL Blazars to the 
Hard X-ray/soft 7-ray Background (shaded area). The 
three SSC curves corresponds to different v pea k values 
(log Vpeak = 12.8, 13.5andl3.8) and are constrained to 
go through the three star symbols representing 1) the to- 
tal contribution of Blazars at 94 GHz (8 x 10~ 6 of the 
CMB intensity), 2) the average 5 GHz to 94 GHz slope 
(a5_94G#2 = 0.2) and 3) the average spectral slope be- 
tween 94 GHz and 1 keV (( a px ) = 1.07). 



4.4. *y-ray Background 

The SSC distributions of Fig. EH predict a negligible 
Blazar contribution to the extragalactic 7-ray Background 
above 100 MeV. Nevertheless, it is well known that Blazars 
are the large majority of the extragalactic 7-ray (E > 
100 MeV) i dentified sources dete cted by the EGRET 
experiment (Hartmanet al.(1999)) aboard the Compton 
Gamma-ray Observatory, and therefore are likely to con- 
tribute in a significant way to the 7-ray Background. 
Indeed Padovani et al.(1993)[ on the basis of a small 
number of sources detected by EGRET, concluded that 
Blazars should make a large fraction, if not the totality, of 
the extragalactic 7-ray background. However, these early 
calculations relied upon a small database and had to as- 
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sume no strong variability, a characteristics that has been 
later demonstrated to be very common in 7-ray detected 
Blazars. 

Instead of considering simple average values of the ra- 
dio to 7-ray flux ratio as in Padovani et al.(1993) we es- 



timate the Blazar contribution to the 7-ray Background 
using the full SED of Blazars scaling it to the integrated 
Blazar flux intensity as calculated with eq. ©. Figure ITT1 
compares the energy distribution of the CMB, CXB and 
CGB to the SED of 3C 279, a well known bright EGRET 
detected Blazar, scaled as explained above. 

Figure ITU shows the large variability of 3C 279 in the 
X-ray and 7-ray band. While at X-ray frequencies the con- 
tribution to the CXB ranges from a few percent to over 
10% in the higher states, the predicted flux at 7-ray fre- 
quencies ranges from about 100% to several times the ob- 
served Cosmic Background intensity. This large excess im- 
plies that either 3C279 is highly non representative of the 
class of Blazars, despite the contribution to the CXB is 
consistent with other estimates, or its duty cycle at 7-ray 
frequencies is very low (see Fig. I12|l . 
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Fig. 11. The CMB, X-ray and 7-ray cosmic backgrounds 
with superimposed the SED of the Blazar 3C279 scaled so 
that its flux at 94 GHz matches the cumulative emission 
of the entire Blazar population (encircled star symbol). 

The same approach can be followed with other Blazars 
detected at 7-ray frequencies. In most EGRET detected 
Blazars the SED of LBL Blazars over-predicts the CGB 
by a large factor. 

A way of quantifying the ratio between the 7-ray in- 
tensity predicted assuming that the source were represen- 
tative of the entire population and the actual Background 
intensity is to use, in analogy with eq. JSJ, a microwave 
(94 GHz) to 7-ray (100 MeV) slope a M7 : 
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Fig. 12. The possible contribution of LBL Blazars to the 
7-ray Cosmic Background (shaded area). Simple SSC en- 
ergy distributions fail to reproduce the observed intensity 
and slope of the 7-ray Background. The observed range 
of 7-ray emission with EGRET, normalized to radio flux, 
indicates that the duty cycle at these energies must be low 
(see text for details). 



We also define an <Vy % = 0.994 as the a M7 of an 
hypothetical source that would produce 100% of the CGB 
if representative of the entire class. This value is based on 
a CGB intensity at 100 MeV shown in Fig. Any real 
source with flatter than ffy-y „ B — 0.994 would 
then predict an integrated flux in excess of the observed 7- 
ray Background if representative of the entire population, 
alternatively its duty cycle must be lower than 100%. 

Tabic 01 lists the main properties of the 34 Blazars de- 
tected by WMAP and by EGRET. Column 1 gives the 
source name, columns 2 and 3 give the Right Ascension 
and Declination, column 4 gives the 5GHz radio flux, col- 
umn 5 gives the 94 GHz microwave flux measured by 
WMAP or estimated from the extrapolation from mea- 
surements in lower WMAP channels or from the litera- 
ture, column 6 gives the minimum and maximum 7-ray 
flux measured by EGRET (excluding lower limits) taken 
from the third EGRET catalog (|Hartman et al.(1999)|), 



11.41 



(7) 



column 7 gives the corresponding a M7 values, column 8 
gives the maximum duty cycle allowed assuming that the 
source is representative of the entire Blazar population, 
columns 9 and 10 give the EGRET and WMAP source 
names. We note that all sources in the list must have a 
duty cycle lower than 100%. 

4.5. TeV Background 

All Blazars so far detected at TeV energies are of the HBL 
type. This is readily interpreted within the SSC scenario 
since only objects where the synchrotron radiation extends 
to near or within the X-ray band can produce a corre- 
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sponding Inverse Compton flux that reaches TeV energies, 
at least assuming a single scattering. 

In the following we estimate the Blazar contribution to 
the TeV background in a graphical way as in the previous 
paragraph for the case of the 7-ray Background but only 
considering the HBL component and not the entire Blazar 
population. From a comparison between the normalization 
of the Blazar radio LogN-LogS and that of the Sedentary 
survey (Giommi et al.(1999)| see also Fig.Owhich is rep- 
resentative of extreme HBLs) , we estimate that these ob- 
jects are about 0.1 % of the Blazar population. 

In Fig. H3 we plot the SED of the well known TeV 
detected Blazar MKN421 normalized at 94 GHz so that 
the flux is scaled to 0.1% of the intensity produced by the 
entire population of Blazars. 

From this figure we see that, despite HBLs are a tiny 
minority, their integrated X-ray flux makes a fair fraction 
of the CXB and that their TeV emission may produce a 
significant amount of extragalactic light, even taking into 
account that the real flux level should be lower than that 
shown in Fig. a t TeV energies since only the objects 
closer than z« 0.2 can be detected at these energies. 

We note however that since extreme HBLs such as 
those of the Sedentary survey are very rare (one object in 
several thousand degrees with flux above a few mJy) the 
extragalactic light at TeV energies should be very patchy, 
associated to single sources, rather than a diffuse light re- 
sulting from the superposition of many unresolved discrete 
sources as in the other cosmic backgrounds. 
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Fig. 13. The CMB, X-ray and 7-ray cosmic backgrounds 
with superimposed the SED of MKN 421 scaled to 1/1000 
of the total Blazar contribution to the CMB since the 
radio LogN-LogS of extreme HBL like this source (in the 
sedentary survey) is about 1/1000 of the LogN-LogS of all 
Blazars. 



5. Summary and discussion 

Blazars are the only type of AGN known to emit non- 
thermal radiation across the entire electromagnetic spec- 
trum, from radio waves to the most energetic 7-rays. In 
some energy bands the power emitted by these sources is 
orders of magnitude larger than that generated through 
the accretion process by other types of AGN. At these 
frequencies, despite their low space density, Blazars are 
the dominant population in the extragalactic sky and con- 
tribute in a significant way to the extragalactic Cosmic 
Backgrounds. 

A deep understanding of the Blazar contribution to 
the Cosmic Background light is becoming an increasing 
necessity as the microwave, the 7-ray and the TeV band 
are about to be intensively explored by a new generation 
of astronomy satellites and ground-based Cherenkov tele- 
scopes. 

The overall Cosmic Background energy has two 
well understood components: the primordial Black Body 
emission peaking at microwave frequencies, or Cosmic 
Microwave Background (CMB) and the X-ray apparently 
diffuse emission arising from the accretion onto super- 
massive black holes in AGN integrated over cosmic time, 
or Cosmic X-ray Background (CXB). We have shown that 
Blazars add a third non-thermal component that at low 
frequencies contaminates the CMB fluctuation spectrum 
and complicates its interpretation, while at the opposite 
end of the electromagnetic spectrum dominates the extra- 
galactic background radiation. 

Our calculations are based on a new, deep Blazar ra- 
dio LogN-LogS that was assembled combining several ra- 
dio and multi-frequency surveys. The integrated radio flux 
from the entire population has been extrapolated to other 
energy bands using observed flux ratios and broad-band 
Spectral Energy Distributions. 

Our results can be summarized as follows. 



results pre- 
findings of 

Giommi & Colafrancesco(2004) The Blazar contribution 



Contribution to the CMB. The 

sented here confirm and extend the 



to the Cosmic Background intensity at microwave fre- 
quencies ranges between 5 x 10~ 5 to 5 x 10~ 6 of the 
CMB, depending on observing wavelength. The conse- 
quences are twofold: i) there is an apparent temperature 
excess of 5-50 /zK, and ii) the CMB fluctuation spectrum 
is affected by a spurious signal that becomes significant 
or dominant at multipole I ss 300 — 600, (see Figs. [7] 
andlHJ- Since Blazars, as all cosmic sources, are proba- 
bly not distributed in a completely random way across 
the sky, the source clustering may significantly increase 
the amount of contamination, particularly at large angu- 
lar scales ( Gonzales-Nuevo et al.(2005)l. 

We also note that the temperature excess causes 
a bias in the statistical distribution of the primordial 
CMB fluctuation spectrum. If not properly removed, this 
non-thermal foreground radiation may greatly complicate 
the detection of primordial non-Gaussianity that may 
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carry important information on the dynamics of the 



inflationary phase (e.g., Peebles(1997) I 



Contribution to the CXB. The contribution to the 
soft X-ray Background has been derived converting the 
integrated radio and microwave Blazar emission into 
X-ray flux using observed flux ratio distributions. Our 
estimated values of 11-12% of the CXB are in very good 
agreement with the independent measurement of 13% for 
radio loud AGN obtained by Galbiati et al.(2004) using 



XMM-Newton data. For the large part this X-ray flux is 
composed of emission due to the end of the synchrotron 
component of HBL objects, while for about one third the 
X-ray flux is due to the flat Inverse Compton spectral 
component of LBL sources. Given the very different X-ray 
spectral slopes of LBL and HBL Blazars, at energies 
higher than 1 keV, the blend between the (steep) syn- 
chrotron HBL component and the (flat) Inverse Compton 
LBL component must change in favour of the latter. Since 
the spectral slope of the Inverse Compton component is 
very similar to that of the CXB up to 40-50 keV, the LBL 
contribution should stay roughly constant around 4-5% 
up to that energy. Above « 50 keV the observed CXB 
steepens, and the contribution of LBL Blazars should 
progressively increase. One of the methods that we have 
used is based on the / x // r distribution of Fig. 0] where 
LBL sources are the large majority among Blazars while 
HBL Blazars are only about 5%. We have converted 
the integrated radio flux into X-rays assuming that the 
/x//r distribution remains constant at all radio fluxes. 
A strong increase of the fraction of HBL Blazars at low 
radio luminosity (hence fluxes) , as required by the Blazar 
sequence within the unified schemes of Fossati et al.(1997) 
and Ghisellini et al.(1998) would result in a much larger 
Blazar contribution to the CXB, inconsistent with the 
observations. This conclusion, together with the find 
ings of |Giommi et al.(2002b)| |Padovani et al. (2003b) 



Caccianiga fc Mar cha(2004)[ Giommi et al.(2005) and 
Nieppola et al.(2005) cast serious doubts on the validity 
of the Blazar sequence. 

Contribution to the soft 7-ray Background. At 

energies higher than «100 keV the contribution to the 
Cosmic Background flux was estimated converting the ex- 
pected Blazar contribution to the diffuse background at 
microwave frequencies using SEDs predicted by simple ho- 
mogeneous Synchro-Self-Compton models, constrained to 
have a microwave-to- X-ray spectral slope a^ x equal to the 
observed average value of 1.07. 

The Hard X-ray/Soft 7-ray Background is reproduced 
both in shape (spectral slope) and intensity assuming an 
average SSC distribution where the synchrotron compo- 
nent peaks at 10 13 5 Hz, well within the range of observed 
values. We conclude that the Cosmic Background between 
«0.5 and « 10 MeV is consistent with being due to the 
tail of the Inverse Compton component of LBL objects 
(see Fig. CHI) . 



Contribution to the 7-ray Background. The SED 

of Blazars detected by EGRET predict much more 7-ray 
Background than observed and therefore they cannot 
be representative of the entire population, at least in a 
stationary situation. Either EGRET detected sources are 
special, non representative objects, or their 7-ray duty 
cycle must be low. Indeed, strong variability at 7-ray 
energies is very common and most objects have been 
detected in widely different intensity states. A plausible 
scenario for the origin of the 7-ray Background is that 
it is due to a mixture of Inverse Compton radiation 
produced by LBL Blazars during strong flares (Fig. I12fl 
and perhaps a less variable component due to the still 
rising part of the Compton spectrum in HBL objects 
(Fig.EJ. 

A TeV cosmic Background? The existence of a 
Cosmic Background at TeV energies has not yet been 
established. Blazars of the HBL type, especially those 
where the synchrotron peak is located at very high 
energies can in principle produce a significant integrated 
flux (Fig. I13fl . A precise prediction is however difficult 
given the uncertainties in the amount of absorption 
of TeV photons via interaction with photons of the 
Infra-red background. In addition, we note that since 
the space density of extreme HBL is very low (less than 
one object in one hundred square degrees at a radio flux 
of 3.5 mJy) the TeV background produced by Blazars 
is composed of widely separated discrete sources rather 
than an apparently diffuse light like e.g. the Cosmic X-ray 
Background. 

In all our estimations we have assumed that the broad 
band SED of Blazars is characterized by approximately 
equal power in the syncrotron and inverse Compton 
components. We note that an average SED, where the 
inverse Compton component is much more powerful than 
the synchrotron component, is not acceptable as it would 
predict a much larger than observed background at hard- 
X-ray/soft Gamma ray energies (see Fig. 1100 . The large 
inverse Compton emission compared to synchrotron flux 
observed in some LBL sources (e.g., Ballo et al.(2002) 1 



should therefore be associated to cither rare objects or 
transient events rather than to an average emission. 

Blazars and future observatories. Figure El shows 
the observed SED of the well known LBL Blazar 3C279, 
scaled down by a factor 1000, overlaid with the limit- 
ing sensitivities of the upcoming Planck and GLAST 
satellites and with a 0.5-10 keV sensitivity that can be 
reached by existing X-ray observatories. This hypothet- 
ical ~ 10 mJy LBL Blazar (or 1 milli-3C279) is at the 
limit of the Planck sensitivity, it is detectable with a 
deep Swift exposure (or a less deep XMM-Newton or 
Chandra observation) and it is detectable by GLAST 
during strong flares. Since the radio LogN-LogS of Fig. 
El predicts a space density of ^ 5 objects per square 
degree with flux above 10 mJy, the Planck mission 
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should detect « 100,000-200,000 Blazars in the « 30, 
square degree high galactic latitude sky. A fractior 
these sources will also be detected by GLAST when f 
ing, exactly how many strongly depends on the duty cy 

A deep all-sky X-ray survey to clean the data fr 
Planck and future CMB missions and to prov 
a database of > 100,000 Blazars for GLAST e 
other 7-ray observatories. Given the significant imp 
of the Blazar foreground emission on the CMB power sj 
trum it is important to remove as much as possible 
contaminating component from the CMB. One poss 
efficient way to achieve this is to exploit the fact that 
spectral slope distribution between microwave and soft 
ray flux of LBL Blazars is very narrow (see Fig. |3J) v 
a dispersion that is probably mostly due to intrinsic v 
ability. The soft X-ray flux of LBL Blazars (that is >£ 
of the Blazar population) is therefore a very good estii 
tor of the flux at microwave frequencies and could be u 
to locate and remove foreground Blazars from the CMB. 
Figures l9l and ITTI imply that a hypothetical all-sky survey 
with limiting sensitivity of a few 10 -15 erg cm _2 s _1 in the 
soft X-ray band, would detect the large majority of Blazars 
above the limiting sensitivity of Planck and therefore allow 
the construction of a database including well over 100,000 
Blazars with flux measurements at radio, microwave and 
X-ray frequencies. Although this type of surveys are cur- 
rently not planned for the near future, they are clearly 
needed to address the previous issues. The extremely large 
sample of Blazars produced by such survey could also be 
used to study with great detail the statistical properties 
of Blazars, including the spatial correlation function and 
would provide a large number of targets for the next gen- 
eration of 7-ray observatories such as AGILE, GLAST and 
future instruments operating in the still poorly explored 
MeV spectral region. 
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Table 4. The list and properties of all WMAP-detected Blazars associated to EGRET 7-ray sources 



Blazar Name 




R.A. 




Dec 




Radio flux 


WMAP flux 


EGRET flux 






Duty 


EGRET 


WMAP 




J2000 


.0 




J2000.0 


5 GHz 


94 GHz 


MOO MeV 






cycle 


name 


catalog 


















Jy 


Jy 


lO^cm-V 1 






% 


3EG J 


number 


(1) 




(2) 




(3) 




(4) 


(5) 


(6) 




(7) 


(8) 


(9) 


(10) 


4C15.05 


02 


04 


50 


3 


15 


14 


10 


3.073 


1.6* 


9-53 


0. 


.846-0.914 


2-12 


0204+1458 


092 


Uy0208-512 


02 


10 


46 


.2 


-51 


01 


02 


3.198 


1.8 


35-134 





.816-0.867 


1-4 


0210-5055 


158 


B2 0234+28 


02 


37 


52 


.3 


28 


48 


08 


2.794 


2.1* 


11-31 


0. 


.877-0.917 


5-13 


0239+2815 


093 


CTA26 


03 


39 


30 


.8 


-01 


46 


35 


3.014 


3.2 


13-178 





.827-0.926 


1-17 


0340-0201 


106 


PKS 0420-01 


04 


23 


15. 


.7 


-01 


20 


32 


4.357 


3.9 


9.3-64.2 





.873-0.946 


4-29 


0422-0102 


110 


Uy0454-463 


04 


55 


50. 


.7 


-46 


15 


59 


1.653 


3.8 


5.5-22.8 


0. 


.911-0.966 


11-47 


0458-4635 


151 


Uy0454-234 


04 


57 


03. 


.1 


-23 


24 


51 


1.863 


2.7 


8.1-14.7 


0. 


.915-0.938 


13-23 


0456-2338 


128 


PKS 0506-61 


05 


06 


44 


.0 


-61 


09 


40 


1.211 


1.1* 


6-29 


0. 


.855-0.915 


3-13 


0512-6150 


154 


Uy0537-441 


05 


38 


51. 


.3 


-44 


05 


11 


4.805 


6.7 


16.5-91.1 


0. 


.880-0.945 


5-28 


0540-4402 


148 


PKS 0735+178 


07 


38 


07. 


.3 


17 


42 


18 


1.812 


1.7* 


15-29 





.872-0.896 


4-8 


0737+1721 
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B2 0827+24 


08 


30 


52. 


.0 


24 


10 


57 
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2.6* 


16-111 


0. 
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0829+2413 
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S50836+710 


08 


41 


24 


.4 


70 


53 


40 
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9-33 


0. 


.854-0.903 
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0845+7049 
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08 


54 


48. 


.8 


20 


06 


30 


2.908 


2.5 


9.7-15.8 


0. 


.910-0.928 
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0853+1941 


115 


4C 29.45 


11 


59 


31. 


.7 


29 


14 


43 


1.461 


2.1 


7.5-163.2 


0. 


.814-0.931 


1-19 


1200+2847 


111 


PKS1221-82" 


12 


24 
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.3 
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13 


10 
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11-36 
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1249-8330 
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12 


29 
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.3 


02 
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04 
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9.0 


8.5-48.3 
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13-73 


1229+0210 
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3C279 


12 


56 


11 


.0 
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47 


19 
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19.0 
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1255-0549 
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PKS 1313-333 


13 


16 


07 
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38 


59 


1.093 
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15-32 





.858-0.887 


3-6 
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12 
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17 


40 
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42 
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PKS 1814-63 c 


18 


19 
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5 
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56 
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.5 


-47 


14 


37 
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1.3* 
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0. 
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3-10 
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BL Lac 


22 


02 


43 


.2 


42 


16 


39 


2.940 


3.8* 


9-40 





.890-0.947 


7-29 


2202+4217 
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PKS2209+236 


22 


12 
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.9 


23 


55 


39 


1.123 


1.3* 


7-46 





.844-0.916 


2-13 
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CTA102 


22 


32 


36 


.3 


11 


43 


50 
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3.1 
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4-18 
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047 
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22 


53 


57. 


.6 


16 


08 


52 


14.468 


5.9 


24.6-116.1 


0. 
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3-16 
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23 


54 
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.6 


45 


53 


03 
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12-43 


0. 
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4-15 
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074 



* Not detected by WMAP in the W-band; flux taken from 90GHz measurements 

** Detected only once by EGRET reported in NED or extrapolated from lower (K,Ka,Q, V-Band) WMAP channels. 

a This source does not have an optical spectroscopic identification yet. However, the radio SED is very flat, typical of Blazars. We 

propose this as a Blazar candidate and associate it with the EGRET source 3EG J 1249-8330 (previously tentatively associated 

to PMNJ1249-8303 by |Mattox et al.(1997)| |. 

6 This EGRET detection has also been associated to the bright radio source PMN J1647-6437 | |Mattox et al.(1997)t which is 
further away from the 7-ray position. 

c This source is a radio galaxy with steep radio spectrum at centimeter wavelengths. At microwave frequencies where WMAP 
is sensitive its spectrum is flat and indistinguishable from a that of a Blazar. 



